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We compute the conserved quantities of the four-dimensional Kerr-Newman-dS (KNdS) black hole 
through the use of the counterterm renormalization method, and obtain a generalized Smarr formula 
for the mass as a function of the entropy, the angular momentum and the electric charge. The first 
law of thermodynamics associated to the cosmological horizon of KNdS is also investigated. Using 
the minimal number of intrinsic boundary counterterms, we consider the quasilocal thermodynamics 
of asymptotic de Sitter Reissner-Nordstrom black hole, and find that the temperature is equal to 
• the product of the surface gravity (divided by 2n) and the Tolman redshift factor. We also perform 

I a quasilocal stability analysis by computing the determinant of Hessian matrix of the energy with 

, respect to its thermodynamic variables in both the canonical and the grand-canonical ensembles 

' and obtain a complete set of phase diagrams. We then turn to the quasilocal thermodynamics of 

four- dimensional Kerr-Newman-de Sitter black hole for virtually all possible values of the mass, the 
rotation and the charge parameters that leave the quasilocal boundary inside the cosmological event 
horizon, and perform a quasilocal stability analysis of KNdS black hole. 
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■ The observational evidence for a positive cosmological constant and the proposed de Sitter conformal field theory 
' (dS/CFT) correspondence ||| have provided motivations for studying the thermodynamics of de Sitter spacetime in 

( <~>) the presence of black holes 3 . Since asymptotic de Sitter black holes have a cosmological event horizon, one may use 
^vq a local description of thermodynamics of the event or cosmological horizons separately . Another way of studying 
, the thermodynamics of these kind of black holes is to investigate their quasilocal thermodynamics P| • 
' One of the quantities associated with the gravitational thermodynamics of black holes is the physical entropy S 
which is proportional to the area of the event horizon(s) The surprising and impressive features of the area law 

^ ^ of entropy is its universality which applies to all kinds of black holes and black strings 0,0, and it also applies to the 
^ [ cosmological event horizon of the asymptotic de Sitter black holes . Another quantity of interest is the temperature 

■ /3"^, which is proportional to the surface gravity of the event horizon(s). 

^ Other black hole properties such as energy, angular momentum and conserved charges, can also be given a ther- 

• ^ , modynamic interpretation [Ifl]. But as is known, these quantities typically diverge for asymptotic flat, AdS, and dS 
' spacetimes. A common approach of evaluating them has been to carry out all computations relative to some reference 

H " spacetime that is regarded as the ground state for the class of spacetimes of interest Unfortunately, it suffers 

■ - - 1 from several drawbacks. The choice of reference spacetime is not always unique p^ . nor is it always possible to embed 
a boundary with a given induced metric into the reference background. Indeed, for rotating spacetimes this latter 
problem creates a serious barrier against calculating the subtraction energy, and calculations have only been performed 
in the slow-rotating regime ^13] . An extension of this approach that addresses these difficulties was developed based 
on the conjectured AdS/CFT correspondence for asymptotic AdS spacetimes 0,0,0], and recently for asymptotic 
de Sitter spacetimes ^3i^l,0|. Indeed, the near-boundary analysis of asymptotically AdS spacetimes can be ana- 
lytically continued to asymptotic dS spacetimes and therefore the counterterms have a straightforward continuation. 
It is believed that appending a counterterm, Ict^ to the action which depends only on the intrinsic geometry of the 
boundary(ies) can remove the divergences. This requirement, along with general covariance, implies that these terms 
are functionals of curvature invariants of the induced metric and have no dependence on the extrinsic curvature of the 
boundary (ies). An algorithmic procedure exists for constructing 1^ for asymptotic AdS '2^ and dS spacetimes T^l, 
and so its determination is unique. Addition of Id will not affect the bulk equations of motion, thereby eliminating 
the need to embed the given geometry in a reference spacetime. Hence thermodynamic and conserved quantities can 
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now be calculated intrinsic for any given asymptotically AdS or dS spacetime. The efficiency of thi s appro ach has 
been demonstrated in a broad range of examples for both the asymptotic AdS and dS spacetimes [13, EM llll I^E 123 ■ 

Recently, we have considered the effects of including I^t for quasilocal gravitational thermodynamics of Kerr, Kerr- 
AdS and Kerr-dS black holes in which the region enclosed by dM was spatially finite. We have also performed a 
quasilocal stability analysis and found phase behavior that was incommensurate with their previous analysis carried 
out at infinity |J,|2J|- this paper, we investigate the thermodynamic properties of the class of Kerr-Newman-dS black 
holes in the context of both temporally infinite and spatially finite boundaries with the boundary action supplemented 
by let- With their lower degree of symmetry relative to spherically symmetric black holes, these spacetimes allow 
for a more detailed study of the consequences of including let for temporally infinite and spatially finite boundaries. 
There are several reasons for considering this. First, one can obtain the thermodynamic quantities corresponding to 
an asymptotically de Sitter spacetime within a fixed radius, while there is no self-consistent thermodynamics for the 
whole spacetime. Furthermore, the inclusion of I^t eliminates the embedding problem from consideration, whether or 
not the temporary infinite limit is taken, and so it is of interest to see what its impact is on quasilocal thermodynamics. 

The outline of our paper is as follows. We review the basic formalism in Sec. m In Sec. IIIII we consider the Kerr- 
Newman-dS4 spacetime and compute the thermodynamic quantities associated to the cosmological horizon. Also we 
give a generalized Smarr formula and investigate the first law of thermodynamics for the cosmological event horizon. 
In Sec. IIVI we consider the quasilocal thermodynamics of Reissner-Nordstrum-dS spacetime, and perform a stability 
analysis through the use of the determinant of Hessian matrix of the energy with respect to the entropy and the 
charge. In Sec. we compute the determinant of Hessian matrix of the energy with respect to its thermodynamic 
variables for Kerr-Newman-de Sitter black holes, numerically, and present our results graphically. We finish our paper 
with some concluding remarks. 



II. GENERAL FORMALISM 

The gravitational action of four-dimensional asymptotically de Sitter spacetimes M, with boundary 5M^ in the 
presence of an electromagnetic field is 

1 f 1 f*-^^ 

loTT 87r Jg^- 

where F^i, — d^Ai, — d^A^ is the electromagnetic tensor field and is the vector potential. In Eq. is the bulk 

manifold with metric g^,y, and SM/^ are spacial boundaries at early and late times with induced metrics "fij and the 
extrinsic curvature O^". The first term is the Einstein-Maxwell term with positive cosmological constant A = 3/P, 
and the second term is the Gibbons-Hawking boundary term that is chosen such that the variational principle is well 

defined. The notation Jg^- indicates an integral over the late-time boundary minus the integral over the early-time 
boundary which are both Euclidean denoted by X+ and T~ respectively. In general, Iq of Eq. ^ is divergent 
when evaluated on the solutions at T^, as is the Hamiltonian and other associated conserved quantities. Rather than 
eliminating these divergences by incorporating a reference term in the spacetime, a new term Id is added to the action 
that is a functional only of boundary curvature invariants [l9l | . Although there may exist a very large number 
of possible invariants, one could add only a finite number of them in a given dimension. Quantities such as energy, 
mass, etc. can then be understood as intrinsically defined for a given spacetime, as opposed to being defined relative 
to some abstract (and non-unique) background, although this latter option is still available. In four dimensions the 
counterterm is |17f 




(2) 

where R is the Ricci scalar of the boundary metric -fij and Jgj^± indicates the sum of the integral over the early and 
late time boundaries. Although other counterterms (of higher mass dimension) may be added to Id, they will make 
no contribution to the evaluation of the action or Hamiltonian due to the rate at which they decrease toward infinity, 
and we shall not consider them in our analysis here. It is worthwhile to mentioning that the inclusion of matter fields 
in the gravitational action requires the addition of further counterterms when the dimension of the bulk is grea ter 
than four. But since we consider the four-dimensional spacetimes, we don't need any additional counterterms |23- 

A thorough discussion of the quasilocal formalism has been given elsewhere [T]| and so we only briefly recapitulate 
it here. The total action can be written as a linear combination of the gravity term Q and the counterterms l(5Jl as 



I^Ig + let- 



(3) 
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Using the Brown and York definition j^J one can construct a divergence- free stress-energy tensor from the total action 
© either on an early or late boundary as p^: 



rp—ab 
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To obtain the boundary stress tensor, we should evaluate Eq. Q at fixed time and send time to infinity. We will 
suppress the -I- sign notation in 7^+"'' in the rest of the paper. To compute the thermodynamic and conserved quantities 
of the spacetime, one should write the metric jab on equal time surfaces in the ADM form: 

-labd-x^dJ" = N'^dp^ + cTy {d(t>' + V'dp) {d(tP + VHp) , 

where the coordinates (jf are the angular variables parameterizing the closed surfaces around an origin. Then the 
total thermodynamic energy of the system is 



Jb 



(6) 



where is the unit normal vector on a surface of fixed p and a is the determinant of the metric aij . If the boundary 
of the spacetime has an isometry generated by a Killing vector then the conserved charge associated with is 



Jb 



(7) 



Since all the asymptotic de Sitter spacetimes have an asymptotic isometry generated by d/dt, there is a notion of a 
conserved mass as computed at future infinity as 



M= f d^(f>y^TabU''^\ 
JB 



(8) 



We compute the quantities of Eqs. (jSJ on a surface of fixed time and then send time to infinity so that it approaches 
the spacetime boundary at . Indeed, for r greater than the radius of the cosmological horizon t is a spatial direction, 
while r will become a temporal coordinate. Thus, the large r surfaces approach I^, and since we purpose to define 
the stress tensor, and the mass at X*, we will evaluate these quantities on the surface of large r. 
Similarly if the surface B has an isometry generated by = d/dcj), then the quantity 



J : 



d^c^V^TabU^-C!" 



(9) 



can be regarded as a conserved angular momentum associated with the surface B. Many authors have used the 
countcrterm method for various asymptotic de Sitter spacetimes and computed the associated conserved charges 
[l7, 18, 19]. We shall also study the implications of including the counterterms ||2J) for Reissner-Nordstrom-dS and 
Kerr-Newman-dS black holes in spatially finite region. 



III. THE THERMODYNAMICS ASSOCIATED TO THE COSMOLOGICAL HORIZON OF THE 

KERR-NEWMAN-DS4 BLACK HOLE 



Here we consider the charged rotating black holes in four dimensions, whose general form is 



A2 / n \2 ^2 



p2 Aesin^e / ,^ (r^ + a^) , , \ 

~d0^ H adt - - — - — -d(j> , (10) 
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where 



Afl = 1 + -7T cos^ 6*, 
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p2 = + 0^008^6', (11) 



and the vector potential is: 



A. = -f.i^l-'^^t)- (12) 

For 9 = 0, the metric H10|l and pi|l is Kerr-dS4 spacetime (or Kerr spacetime if I oo), and for a = the metric 
is Reissner Nordstrom-dS4 spacetime, which has zero angular momentum. Here the parameters m, a, and q are 
associated with the mass, the angular momentum and the charge of the spacetime, respectively. The metric of Eqs. 
(|10|l and Hll|l has three horizons located at r_, r_|_, and Tc, provided the parameters to, Z, a, and q are chosen suitably. 
The three parameters I, a and q should be chosen such that l'^ + — lAa^P — 12q^P > 0, and m should lie in the 
range between TOj^^j-j^ < to < TO2gj-it, where m^f,j.^^ and TOj^j-jt are the positive real solutions of the following equation: 

27^^771^ + (33Z^a^ - 36l^q^ - SSZ^a^ + 361'^q^a^ + Pa^ - f)rn? 

+l2aH^q^ + 4/2a8 + + + + 8l\^ + 8l\^a^ = 0. (13) 

It is worthwhile to mention that when m = m^^j.^^, r_ and r+ are equal and we have an extreme black hole. For the 
case in which m = 'rijcriti the radii of event and cosmological horizons are equal and we have a critical hole. Also, 
one may note that in the limit I —> oo, TOj^^j-j^ = \fcP^-\-cp' ^ ™2crit ~^ J^i = to ± \J m? — — q^ , and Tc — > oo. 

For given values of a and to, the parameter I should be greater than ^cj-iti where If^^-n is the positive real solution 
of Eq. (fT!^ . When / = Zgj.j^, we have a critical hole. Also the positive real solution of Eq. (fT!^ for a gives the critical 
value of a (upper bound for a) for which the hole is extremal. The two real solutions of Eq. H13() for q determine the 
lower and upper bounds of q that belong to the critical and extremal black hole, respectively. 

Since the temperature associated with the event or cosmological horizons has different values for asymptotic de 
Sitter spacetimes one should use a local description of the thermodynamics of the black hole. In the remain of 
this section we investigate the thermodynamics associated to the cosmological horizon and in the next section we 
investigate the thermodynamic quantities associated with the spacetime within an arbitrary surface with fixed radius 
r. 

Analytical continuation of the Lorentzian metric hy t —t it and a ~* ia yields the Euclidean section, whose regularity 
at r = Tc requires that we should identify (r, </)) ~ (t + /3c, + iPc^c), where Bn and flc are the inverse Hawking 
temperature and the angular velocity of the cosmological event horizon given as j9( : 

P = ^-l'r^irl + -') (14) 

Sr^ + a^r^^-Pr^^+a-^P + qH^' ^ ' 

nc = (15) 

Since the area law of the entropy is universal, and it can be used for the cosmological horizon of dS black holes |3, 
the entropy is 

S^^=-J!i±^, (16) 

As we discussed in Sec. m one can compute the finite action through the use of boundary counterterms. To do 
this one should take a space-like hypersurface, V of large constant r (r^ < r < oo) in the future and calculate the 
total action In the limit that r tends to infinity, this surface tends to future infinity X"*" as expected. The first 
integral of the action in ^ should be first computed from rc to r, and then the total action should be integrated on 
the hypersurface V. Sending r to infinity the total action of the system is : 

, ^ (rg + a^)(rg + r^a^ + mP) - r^Pq^ 
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where and j3c are the radius and the inverse of the Hawking temperature of the cosmological event horizon. The 
total mass can be calculated through the use of Eq. © and sending the temporal coordinate to infinity. We obtain 

M = -g. (18) 

It is worthwhile to note that as a and q go to zero, the action and mass given by Eqs. (|17|l and (|18|l reduce to those 
of Schwarzschild-dS black hole obtained in |l7j |. The total angular momentum can be computed through the use of 
Eq. ^ for a boundary with finite radius, r+ < r < rc- It is a matter of calculation to show that J is independent of 
r given as: 

ma 

J = -^. (19) 

We now obtain the mass as a function of the extensive quantities S, J and Q, where Q defined as (//S. Using the 
expressions ((T51l and (fTB|l for the mass, the angular momentum and the entropy, and the fact that Ar(rc) = 0, one 
can obtain a generalized Smarr formula as 

It is worthwhile to mentioning that the analytic continuation / — + il (Wick rotates) changes the above Smarr formula 
to the Smarr formula of the Kerr-Newman-AdS black hole given in Ref. Now, one may regard the parameters 
S, J and Q as a complete set of extensive parameters for the mass M{S,J,Q) and define the intensive parameters 
conjugate to S, J and Q. These quantities are the temperature, the angular velocities, and the electric potential 
defined as 

It is a matter of straightforward calculation to show that the temperature calculated from Eq. (|21|l coincides with 
Eq. H14|) . and the angular velocity obtained from Eq. I|21l) is the thermodynamic angular velocity that is equal to 
ilc ~ a/P. Also one can obtain the electric potential <I> as 

<I' = -^^. (22) 

It is worthwhile to note that the thermodynamic quantities calculated in this section satisfy the first law of thermo- 
dynamics, 

dM = TdS + ndJ + ^dQ. (23) 



IV. THE QUASILOCAL THERMODYNAMICS OF A REISSNER-NORDSTROM-DS4 BLACK HOLE 

First we investigate the quasilocal thermodynamics of a non-rotating case in which all the thermodynamic quantities 
can be integrated easily. For a = 0, one obtains the quasilocal energy of a Reissner-Nordstrom-dS black hole as 

If r^X I 2m , , 

^ = 2('-Vj-V^-]^-v + ?^- (^4) 

Using the expressions Q = q for the charge and S = 7r(r^ + for the entropy of the event horizon, one can 

write E as a, function of thermodynamic quantities S and Q. Now identifying the energy (|24|l as the thermodynamic 
internal energy for the Reissner-Nordstrom-dS black hole within the boundary r, then the corresponding temperature 
and electrostatic potential at the boundary with fixed radius rj^<r<rc are 
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FIG. 1: [d'^E/dS'^)Q versus r+ iov I = 1, a = 0, q = 0.2, (solid), and (dotted). 
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FIG. 2; Determinant of Hessian matrix (solid) and {d'^ E /dS^)Q (dashed) versus m for I — 1, a = and q — 0.2. 



The first factor in the above expressions is the inverse of the lapse function for Reissner-Nordstrom-dS4 metric 
evaluated at r, and is the Tolman redshift factor for temperature in a static gravitational field 2^. The second 
factors in Eq. H25|) and (|25|l are the surface gravity (kh) of the black hole divided by 27r, and the electrostatic 
potential of the hole relative to electrostatic potential at r. Therefore as in the case of Schwarzschild-AdS and 
Reissner-Nordstrom-AdS black holes [28L the temperature corresponding to the spacetime within the surface of 
fixed r is the product of kh/^tt and the redshift factor. Note that the term due to the counterterm does not affect 
the temperature and the electrostatic potential of the black hole. 



A. Stability analysis in the canonical and the grand canonical ensemble 

The local stability analysis in any ensemble can, in principle, be carried out by finding the determinant of the Hessian 
matrix of S with respect to its thermodynamic variables, H^.^^^ — \d'^S/dXidXj\ where AT^'s are the thermodynamic 
variables of the system. Indeed, the system is locally stable if the determinant of the Hessian matrix satisfies H^. ^. < 
[30l |. Also, one can perform the stability analysis through the use of the determinant of Hessian matrix of the energy 
with respect to its thermodynamic variables, and the stability requirement x — ^ may be rephrased as Hy y > 

El- 

The number of the thermodynamic variables depends on the ensemble which is used. In the canonical ensemble, 
the charge is a fixed parameter, and therefore the positivity of the thermal capacity {d'^E/dS'^)Q is sufficient to assure 
the local stability. It is a matter of calculation to show that {d'^E/d^S)Q is 

~2r+(3r4 - + rlP){Pr^ - r* + q^l% (27) 

where Aq is at a = 0. Note that Aq is equal to zero at horizons, and therefore {d^E/d^S)Q goes to infinity as 
r approaches the radius of event or cosmological horizon, we set the radius of the boundary equal to r — 0.5. For 
a given I and a fixed value of r in the range r+ < r < r^, Eq. H27fl has two real solutions for r+, which shows that 
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there exist two stable phases separated by an unstable intermediate mass phase. This is not true for the case of an 
uncharged hole for which Eq. 127|l has a single real solution (see Fig. 

In the grand-canonical ensemble, the stability analysis can be carried out by calculating the determinant of Hessian 
matrix of the energy with respect to S and Q, 

l6TT^r-^{r - r+) ^ + + + 

+2r+(f - r^){rPq^ - 2r+Pq^ - rl^l - Srr^ )} 

x{rr+^ + r'^r+^ - r+fr + r+r^ + q^l^Y"^. (28) 

Again note that the determinant of Hessian matrix goes to infinity as r goes to the radius of the event or cosmological 
horizon. Since the number of thermodynamic variables in the canonical ensemble is more than that of the grand- 
canonical ensemble, the region of stability is smaller for the latter case. This fact can be seen in Fig. |2 which displays 
{d^E/d^S)Q and the determinant of the Hessian matrix as a function of the the mass parameter m. 



V. THE STABILITY ANALYSIS OF THE KERR-NEWMAN-DS4 BLACK HOLE 

Now we turn to the quasilocal thermodynamics of Kerr-Newman-dS4 black hole. The energy of the system can be 
calculated through the use of Eq. ((BJ where the boundary S is a two dimensional surface with fixed radial coordinate 
r+ < r < Vc- This energy can be identify as the thermodynamic internal energy for the Kerr-Newman-dS spacetime 
within the boundary r. 

It is important to say that the total angular momentum calculated by Eq. Q does not depend on the radius of the 
boundary surface and it is equal to the value given in Eq. (|19l) . Now using Eq. H19() for the angular momentum, and 
the expressions Q = qfE. and S = T^{r'\_ -t- a^)/4 for the charge and the entropy, one can obtain the metric parameters 
in terms of S*, J and Q: 

^ = 2Z2^3/2gi/2r (l + 4^'^'^'^r'')'' (29) 

where 

r = 54 - An^J^l^S + A^\JH^ + tt^QH* + 27r3g2;45 _ 27r2Q2/2^2 ^ ^2^4^2 _ ^S'^^. (30) 

With these expressions for a, m and q, we can write as a function of thermodynamic quantities S, J and Q, and 
regard E{S, J, Q) as the thermodynamic internal energy for the Kerr-Newman-dS black hole within the boundary 
with radius r. In the canonical ensemble, the charge and angular momentum are fixed parameters, and for this 
reason the positivity of the {d'^E/dS'^)j^Q is sufficient to assure the local stability. For general values of the metric 
parameters, however, we cannot analytically calculate the internal energy as a function of S, J and Q, and so we 
resort to numerical integration. To do this, we find e in Eq. @as a function of S, J and Q, calculate the derivative 
of e{S, J, Q) and then compute the integration numerically. We set the radius of the boundary equal to r = 0.56, and 
we apply the restrictions discussed in Sec. Illll on the parameters I, m, a and q. Note that for r = 0.56, each of these 
parameters has two critical values, one due to the extreme Kerr-Newman-dS black hole and the other due to the case 
in which this radius is equal to the radii of event and cosmological horizons. 

Figure |31 displays the m dependence of {d^E/dS^)j^Q and shows that for a given value of a and a large value of 
q, only a stable phase exists for all the allowed values of m (wicrit < ™ < '^2crit)' 1 decreases there begins 

to appear an unstable phase between the two stable phases. Also Fig. 0]shows the m dependence of [d^E / dS^)j^Q 
for small q and different values of a. It points out that there exists an unstable phase, between two stable phases, 
which begins to disappear as a increases. Figures |5l and |H1 display the a dependence of {d^E/dS^)j,Q. They show 
that for large values of m and q the hole is stable while for large m and small q, there exist two stable phases 
separated by an unstable intermediate angular momentum phase. Also note that for small and intermediate q or 
m, the low angular momentum phase is not stable (see the dotted curve in Fig. |31and|njl. Figure [7| which displays 
the q dependence of {d^E/dS^)j^Q shows that for a large value of m and a there exists only a stable phase, but an 
unstable phase begin to appear as a decreases. We also find that for small a and intermediate m the black holes with 
small q are not stable (see Fig|Hl. It is worthwhile to note that in all of these cases {d^E/dS^)j^Q goes to infinity as 
the hole approaches its extreme case and therefore the nearly extreme black hole is stable in the canonical ensemble. 



8 



4- 












3- 












2- 
1- 















0.12 0.14 a.16 0;T8 0.2 


0.22 



FIG. 3: {d'^E/dS^)j,Q versus m for / = 1, a = 0.1, q = (dotted), 0.1 (thick-line), and 0.2 (solid). 
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FIG. 4: {d^E/dS^)j,Q versus m for / = 1, g = 0.1, a = 0.1 (thick-line), and 0.15 (sohd). 

This stability analysis is in qualitative agreement with that of Davies 32], who performed a thorough analysis of 
the thermodynamic properties of asymptotic Kerr-Newman-dS black holes. However, Davies used the ADM mass 
parameters, whereas we are considering the quasilocal energy E in Eq. lO as a function of r, 5, J and Q. 

In the grand-canonical ensemble, the thermodynamic variables are the entropy, the charge and the angular 
momentum. The determinant of the Hessian matrix can be calculated numerically. Some typical results of these 
calculations are displayed in Figs. I9I12I As one can see from Fig. I^which shows the m dependence of the determinant 
of Hessian matrix, for a given value of a and large values of q, the hole is stable, but for small values of q, only the 
extreme case is stable. The q dependence of the determinant of the Hessian matrix are shown in Figs. IIUI and ITTI 
Figure [TUI shows that for large values of m and a we encounter a single stable phase, while for large m and small a the 
only stable phase is the nearly extremal hole. We also find that for small a the hole is stable for large m and unstable 
for small m. For intermediate m only the nearly extreme hole is stable (see Fig. Ill|l . The a dependence of the 
determinant of the Hessian matrix (Fig. I12|l shows that for large values of to, when q is large the only stable phase 
is the nearly extreme hole, but as q decreases a stable low angular momentum phase will appear. The comparison 




FIG. 5: {d^E/dS'^)j,Q versus a for Z = 1, m = 0.2, q = 0.2 (thick-hne), 0.15 (dotted), and 0.1 (sohd). 
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FIG. 6: {d'^E/dS^)j,Q versus a ioi I = I, q = 0.1, m = 0.15 (dotted), and 0.2 (solid). 



4- 








3- 








2- 
1- 











0.02 0.04 S.WO:08 'T)T U.T2Tj.14 0.16 0.18 



FIG. 7: {d'^E/dS^)j,Q versus g for Z = 1, m = 0.2, o = 0.2 (solid), and 0.1 (dotted). 
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FIG. 8: {d^E/dS'^)j,Q versus g for / = 1, a = 0.1, m = 0.1 (solid), 0.15 (thick-line), and 0.2 (dotted). 




FIG. 9: Determinant of Hessian matrix versus m for Z = 1, a = 0.1, q = 0.1 (thick-line), and 0.2 (dashed). 
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FIG. 10: Determinant of Hessian matrix versus q for I — 1, m = 0.2, a = 0.2 (dotted), and 0.1 (solid). 
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FIG. 11: Determinant of Hessian matrix versus q for I — 1, a — 0.1, m = 0.15 (dashed), 0.2 (thick-line), and 0.22 (solid). 



of the figures of the canonical and grand-canonical ensemble shows that the region of stability is smaller in the 
latter case as stated before, and there are cases that even the extreme hole is not stable in the grand-canonical ensemble. 



VI. CLOSING REMARKS 



In this paper, we used the counterterm method to calculate the conserved quantities and the Euclidean action of 
the four-dimensional Kerr-Newman-dS black hole. Also we obtained a generalized Smarr formula for the mass as a 
function of the entropy, the angular momentum and the electric charge and showed that these quantities satisfy the 
first law of thermodynamics associated to the cosmological horizon. Since there is no well-defined thermodynamics for 
the asymptotically de Sitter black holes, a local investigation of the thermodynamics is necessary. Our investigation of 
the boundary-induced counterterm prescription Q at finite distances has allowed us to obtain a number of interesting 




FIG. 12: Determinant of Hessian matrix versus a for I — 1, m = 0.2, q — 0.15 (dotted), and 0.1 (solid). 
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results, which we recapitulate here. 

We were able to compute the quasilocal energy, the temperature, and the Hessian matrix of the energy with respect 
to S and Q for arbitrary values of the parameters of the Reissner-Nordstom-dS black hole analytically. We found 
that the temperature at fixed r is the product of Kn/'in and the redshift factor. We perform a stability analysis both 
in the canonical and grand-canonical ensemble, and obtained a complete phase diagrams. We found that at a fixed 
value of r there exists two stable phases separated by an unstable phase with intermediate mass value. That is, charge 
causes the low mass black hole to be stable, a fact which is not true for the Schwarzschild-dS black holes. 

We also computed the determinant of the Hessian matrix, quasilocally, for arbitrary values of the parameters of the 
Kerr-Newman-dS4 solution, apart from the mild reality restrictions considered in Sec. IIIII These quasilocal quantities 
are intrinsically calculable numerically at fixed r without any reference to a background spacetime. By finding the 
parameters a, m and q in terms of the thermodynamic quantities S, J and Q, we regarded the interior of the quasilocal 
surface as a thermodynamic system whose energy E is a, function of S, J and Q. We studied the stability of the Kerr- 
Newman-dS black hole both in the canonical and grand-canonical ensembles. Thereby we encountered an interesting 
phase structure. We found that in the canonical ensemble the extreme black hole is stable for all the allowed values of 
the metric parameters while this is not true in the grand-canonical ensemble. The stability analysis in the canonical 
ensemble yielded results that are in qualitative agreement with previous investigations carried out by using the ADM 
mass parameter [s^- This stability analysis is, of course local, and phase transitions to other spacetimes of lower free 
energy might exist. 
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